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PREFACE 

This  short  summary  or  oceanic  internal  waves  and  their  effect  on  the  speed  of  sound  was 
prepared  as  an  appendix  to  an  inteiim  report  concerned  with  the  influence  of  internal  waves  on 
the  performance  of  acoustic  surveillance  systems.  Because  of  the  limited  distribution  of  the 
report  it  seems  worthwhile  to  make  the  material  available  separately.  Hopefully,  it  will  assist 
the  reader  in  understanding  the  many  papers  which  have  recently  appeared  dealing  with  the 
subject  of  acoustic  fluctuations  and  ocean  variability.  It  is  concerned  with  deriving  as  rapidly 
and  as  directly  as  possible  the  mathematical  expression  for  the  correlation  function  of  the 
sound-speed  fluctuations  resulting  from  internal  waves.  As  a  regrettable  consequence,  a  great 
deal  of  the  theory  of  internal  waves  is  left  by  the  wayside. 


Din',  ribution. 
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SOME  BASIC  ASPECTS  OF  OCEANIC  INTERNAL  WAVES  AND 
THE  SOUND-SPEED  FLUCTUATIONS  THEY  GENERATE 

Internal  waves  are  those  oscillations  in  the  ocean  which  result  from  a  stable  density 
(entropy)  stratification.  They  are  thought  to  be  responsible  for  motions  with  horizontal  scales 
ranging  from  kilometers  to  tens  of  kilometers,  vertical  scales  ranging  from  tens  of  meters  to 
hundreds  of  meters,  and  temporal  scales  ranging  from  minutes  to  hours.  Experimental  studies' 

7  i  •  V 

conducted  over  a  period  of  years  indicate  internal  waves  exist  in  all  the  world's  oceans,  have  a 
characteristic  (universal)  spectrum  which  is  stationary,  horizontally  homogeneous,  and,  to  a 
large  degree,  horizontally  isotropic. 

Two  developments  have  recently  occurred  which  greatly  enhance  our  ability  to  discuss,  in 
quantitative  terms,  the  relationship  between  internal  waves  and  acoustic  fluctuations.  These 
important  advances  are  (1)  the  construction  in  1972  of  a  general,  phenomenological  inter¬ 
nal-wave  model  by  Garrett  and  Munk2  and  (2)  the  establishment  in  1975  by  Munk  and 
Zachpriasen3  of  the  precise  relationship  between  the  quantities  which  characterize  the  internal- 
wave  field  and  the  induced  change  in  the  sound-speed  profile.  In  this  report  these  two  develop¬ 
ments  will  be  reviewed.  To  this  aim  we  only  consider  those  aspects  of  the  Garrett-Munk  (GM) 

*  »  *  ^  .  *  *  *  .  •  ‘ 
model  which  are  relevant  for  acoustic  signal  coherence  studies. 


We  begin  with,  the  linear  equations  which  govern  internal  waves,  i.e., 
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— (0,Uy  4- 

dt  y 


Po 


Bp,w 

dx 


~  u,u,  +  — 

Po 


Bpiw 

By 


Manuscript  submitted  August  24,  1979. 


0 

0, 


(la) 

(lb) 


1 


D.  R  PALMER 


where 


£  +  iTr  +  ,’w‘-*- 


dux  9uy  fly 

8x  dy  8z 


0, 


(lc) 

(ld) 

(le) 


ux,uy  “  horizontal  components  of  particle  velocity, 
v  -  vertical  component  of  particle  velocity, 

Pi w  —  internal-wave  induced  pressure, 

£  —  vertical  particle  displacement, 

p0  —  equilibrium  density  (taken  to  be  constant), 
n(z)  —  bouyancy  or  Brunt-Vaisala  (B-V)  frequency, 
at,  —  inertial  frequency, 

—  (2  cycles/day)  sin  (latitude), 

—  0.042  cph  at  30°  latitude. 

The  coordinate  system  is  orientated  so  that  the  z  axis  is  positive  downward,  the  x  axis  (the 
direction  of  ux)  is  to  the  east  and  the  y  axis  (the  direction  of  uy)  is  to  the  south. 

These  equations  are  the  result  of  several  approximations.  One  starts  with  the  source-free 
hydrodynamic  equations  valid  for  a  rotating  earth.  The  equations  are  then  linearized  and  the 
Boussinesq  approximation  is  made.  It  is  then  assumed  that  the  cosine  of  the  latitude  may  be 
set  equal  to  zero.  (This  last  approximation  is  called  the  "traditional'  approximation.)  The 
details  of  the  derivation  of  Eqs.  (1)  may  be  culled  from  the  book  by  Phillips4  or  the  one  by 
Eckart.5 
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Equations  (1)  are  linear,  coupled  differential  equations  and  can  therefore  be  solved  by 
Fourier  transforming  in  time  and  then  using  the  technique  of  separation  of  variables.  The  gen¬ 
eral  solution  is 


ux(x,  /) 

ux 

Uy(%,  t) 

-  2n 

uy 

v(x,  r) 

C(x,  t) 

“  Z  J  J  JJrrT  exPl'kjrcos(>l>  -  <t>)  -  iut) 

V 

D 

P/wCx.,  t)/p0 

P 

where 


(2a) 


ux 

1  cwy 

—  Zj(z)(i cos  0  +  —  sin  &) 

Kj  co 

1  Ol, 

—  Zv(z)(/sin  $ - cos 

kj  o> 

Uy 

V 

* 

Zj(z) 

~t ZjU) 

D 

/(<u2-a>,2) 

z'l!) 

P 

Here  x  —  rcos  <t>,  y  —  rsin  <t>  and  the  Aj’s  are  constant  coefficients  (i.e.,  they  are  independent 
of  x  and  t).  The  normal-mode  eigenfunctions  Zy(z)  and  corresponding  eigenvalues  Hr,  are 
determined  by  the  eigenvalue  equation 


Z/(z)  +  k} 

Here,  and  in  Eq.  (2b),  a  prime  indicates  differentiation  with  respect  to  z.  The  eigenfunctions 
are  assumed  to  vanish  at  the  surface  (i.e.,  the  top  of  the  main  thermocline)  and  at  the  bottom 
of  the  ocean.  The  index  J  is  discrete  and  is  taken  to  run  from  one  to  infinity. 

For  sound  coherence  studies  the  effect  of  the  vertical  displacement  is  more  important 
than  that  of  the  pressure  or  the  components  of  the  velocity.3  6  Consequently,  we  need  only 
consider  the  expression  for  { 


n2(z)  —  <o2 


<u2  —  <u2 


Zj(z)  -  0 


(3) 
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£(x,  t)  -  T  f  -~^r  f  77^-  I—  <4y(</»,  at)  exp[/*,/xos(0  -</>)-  /to/]  Zy(z).  (4) 

y  (2w)  -J  (2tt)  (  <0 

In  the  GM  model,  the  amplitudes  ,4,(i/i,a>)  are  taken  to  be  zero-mean,  Gaussian  random  vari¬ 
ables.  Hence,  all  moments  of  the  A/ s  can  be  expressed  in  terms  of  the  correlation  function 
< w)Aj  *  (■/»',  o»')>.  (It  is  possible  to  use  the  reality  of  {  to  relate  <AjAr>  to 
<AjAj  *>)  The  requirements  of  reality,  stationarity,  and  horizontal  homogeneity  dictate  that 
this  correlation  function  have  the  general  form 

<Aj(\li,  <o)Aj*  (ij/', “  (2tt)28(o)  —  oj')8(t|/  —  j(<l>,  to)  (5) 

where  T,  is  a  real  function  having  the  property 

r y(0,  to)  —  Ty(l/»  ±  n,  —to). 

If  the  additional  assumption  of  horizontal  isotropy  is  made,  ry(i/»,  at)  is  independent  of  the 
angle  <// 

r y(l//,  to)  —  r y  (to)  “  Ty(— to)  .  (6) 

With  Eqs.  (4)-(6)  we  have 

oo 

<£(x,  t)C(x',  /')>-  —  L  /  - T  r,<to)y0(*,/Ocos[to(r  -  /')]  Zj(z)Zj(z')  .  (7) 

IT  j  0  (0 

Here  R  -  V(x  -  x')2  +  (y  -  y')2  and  J0  is  the  common  Bessel  function.  We  have  used  the 
fact  that  ry,  Zy,  and  kj  are  all  even  functions  of  to  to  write  the  integral  as  one  over  positive  fre¬ 
quencies  only.  Since  the  A/s  are  Gaussian  random  variables,  so  is  {,  and  all  moments  of  {  can 
be  expressed  in  terms  of  the  correlation  function  given  in  Eq.  (7). 

Garrett  and  Munk  assumed  a  "canonical  form"  for  the  buoyancy  frequency, 

n(z)  -  n0e-,IB.  (8) 

which  is  thought  to  give  a  reasonable  representation  below  the  mixed  layer.  Typically,  n0  “  3 

cph  and  B  -  1.3  km.  They  then  solved  for  the  modal  eigenfunctions  and  eigenvalues  usir.?  the 

WKB  approximation  with  the  result  that 
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z'w  -  V& 


sin  Urnf'z—  ;  u,  <  a>  <  n(z), 
»o 


=  0;  otherwise. 


where  the  C.’s  are  normalization  constants  and 


kj  =*  (w2  -  <u,2)l/2.  (10) 

ano  ,  \ 

By  assuming  that  |z  -  z'\/B  is  small  and  that  short-period  oscillations  associated  with  variations 
in  z  -  y  (z  +  z')  are  unobservable  and  can,  therefore,  be  averaged  out,  one  has 


Zj  (z)  Zy(z’)  =  — ■  cos  [kj  (z)  (z  -  z')),  (11) 

n(z) 

for  w,  <  u  <  n(z)  and  zero  otherwise.  The  quantity  is  called  the  local  vertical  wave 
number  and  is  given  by  the  expression 


*y<2>  ~ 


Substituting  Eq.  (11)  into  Eq.  (7)  gives 


<{(x,  t)  (  (x‘,  /')>  - 


- (zr  Z  fU>  CJ  FjMJoikjR)  cos  M/-r')]  cos  [*y(z)(z  -  z')].  (13) 

n(z)  j  ‘  id 

The  mean-square  vertical  displacement  is  obtained  from  Eq.  (13)  by  setting  x  -  x'  and 
i  -  i.e., 

<({(*,  0)2>  -  -^r  ZS"  CjTjio,).  (14) 

n(z)  y  <u2  J  1 

Because  the  internal-wave  field  is  stationary  and  horizontally  homogeneous  <({(x,  r) ) 2 >  does 
not  depend  on  t  or  on  the  horizontal  coordinates  x  and  y.  It  does  depend  on  z,  however, 

reflecting  the  vertical  inhomogeneity  of  the  internal-wave  field.  By  setting 


<({(*,  f))2>  s  <{2(z)> 
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we  can  simplify  the  notation  somewhat.  It  turns  out  that  C,  r,(ca)/to2  rapidly  decreases  with 
increasing  u>  so  it  is  usually  a  valid  approximation  to  set  the  upper  limit  of  the  integral  in  Eq. 
(14)  equal  to  infinity.  Hence 

<{■(,)>  -  fe’fy1..  os) 

n(z) 

i.e.,  the  product  n(z)<{2(z)>  is  approximately  independent  of  depth.  This  equation  can  be 
used  to  define  an  extrapolated  mean-square  displacement  <(o  >. 


n(z)  <{2(z)  >  -  n0  <£o  >• 

By  fitting  data  Garrett  and  Munk  found 

<{<?  >l/2  -  7.3  m. 

Returning  to  Eq.  (13),  if  we  define  a  normalized  spectrum  by  the  relation 

(Cj/trn (z))w~2  T Aw) 

F,(to,  n  (z))=  - r— - ;  w,  <  <o  <  n(z), 

J  <iHz)> 

—  0;  otherwise 


we  have,  with  z  —  y(z  +  z'), 


(16) 


(17) 


(18) 


<£(x,  r)£(x',  t')>  - 

<£2(z)  >  £  fQ  duFj(u),  n(z))J0(kjR)  cos  Ikj(z)  (z  -  z')]  cos  Mr  -  r')],  (19) 

j 

and,  clearly, 

ltJ~du,FJ('o.n(z))-\  .  (20) 

j 

The  vertical  displacement  spectrum  Fl  (a>,  j\  z)  defined  by  Munk  and  Zachariasen3  is  related  io 
our  F/iw,  n(z))  through 


f((a),7';z)  -  <{2(z)>  Fj  (to,  n(z)). 


(21) 


Consequently 


<{(x,  r){(x‘,  /’) >  - 

*n(F) 

£  Jw  dw  F((w,  j\  z)J0(kj  R)  cos  U,(z)(z  —  z')]  cos  lw(f  -  /’)), 


(22) 
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<C  !W>-Xf  dcF^w.j^z).  (23) 

,  J  ' 

The  integration  in  Eqs.  (22)  and  (23)  extends  from  w,  to  n(z)  as  a  result  of  the  definition  Eq. 
(18). 

In  GM-type  models,  Fj  is  the  product  of  two  factors 

Fj  (<o.n)  -  H(J)G(<o,n)  .  (24) 

(We  have  explicitly  indicated  the  dependence  of  G  on  n  in  writing  Eq.  (24).  In  most  papers 
this  dependence  is  suppressed.)  It  follows  from  the  definition  of  Fj  that 


I  mj)  - 1 , 

i- 1 


f  diD  Giai.n)  —  1  . 
•'0 


The  function  G  is  taken  to  have  the  form 


CD,  (co^-'w/) 

G(omi)  -  Nc  ;  (Oi<w<n  , 

(O 

r 

-  0  ;  otherwise. 

The  normalization  constant  Nc  is  determined  from  Eq.  (26), 


*c--  n  +  o  —  = 


«  IT 


There  are  (at  least)  three  versions  of  the  GM  model,  GM72,2  GM75,7  and  GM75y.‘  They 
differ  in  that  they  assume  different  forms  for  H(J).  The  latest  version,  GM75y,  assumes 


hw  -  ».  -fa  • 
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where9 


Typically,  j  *  -3. 


Nii' 


7-1 


j2+jl 


_1 

27 


7  (w-y.-D  . 


(30) 


Summarizing 


<f(x,  r)£(x'.  r')>  - 


„(F>  <t0> 


X  £  J  dot  H(j)G(<i>,n(z))J0(kjR)cos[Kj(.z)(.z-z')]cos[a>(t-t')\ 

i- 1  0 


where 


z  -  y  (z+z'),  /?  -  V(x-x')2+(.y-.y')2 
n(z)  -  n0e~z/B  , 
k,  -  (<^2-‘-12)1/2  . 

noU 

*y(*>  ”  • 

,  x  4  (to2— to,2)1^2 

G(to,n)  —  —  to,  - : - ;  to,<to<n  , 

it  or 

-  0  ;  otherwise  , 


and 


Typical  values  are 


H(J)  - 


(y2+yi) 


I  02+7  2.)-' 

7-1 


(31) 

(32a) 

(32b) 

(32c) 

(32d) 


(32e) 

(320 


to,  —  0.042  cph  , 
n0  ”  3  cph  , 

B  —  1.3  km  , 

<4o>,/2  ”  2.3  m  , 

J*  -3  . 

Before  proceeding  we  would  like  to  mention  a  caveat.  Much  of  the  theory  of  internal 
waves  outlined  here  is  dependent  on  the  assumption  that  the  bouyancy  frequency  is  a  mono- 
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tone  decreasing  function  of  depth  below  the  mixed  layer.  There  are  situations,  however,  where 
the  bouyancy  frequency  has  a  pronounced  minimum  which  tends  to  trap  the  internal  waves  in  a 
duct.  This  ducted  propagation  can  have  a  pronounced  effect  on  the  character  of  the  acoustic 
signal  fluctuations.10  For  such  a  situation  Eq  (9)  is  not  correct  and  one  might  need  to  actually 
calculate  the  internal-wave  normal-mode  eigenfunctions  numerically"  in  order  to  correctly 
describe  the  sound-speed  fluctuations. 

We  now  wish  to  record  the  relationship  between  the  vertical  particle  displacement  {  and 
the  change  in  the  sound  speed.  In  the  absence  of  internal-wave  activity  the  sound  speed  is 
assumed  to  be  a  non-random  function  of  depth  alone 

c(x,  t)  “  <c(x,  /)>  -  c0(z).  (33) 

When  internal  waves  are  present  the  sound  speed  acquires  a  small  fluctuation  8c  which  is  a  ran¬ 
dom  function  of  all  three  spatial  coordinates  and  of  the  time 

c(x,  i)  -  c0(z)  +  8 c(x,  r).  (34) 

By  definition  Sc  has  a  mean  value  of  zero 

<8c(x,  r)  >  =  0. 

Munk  and  Zachariasen3  argued  that  Sc  is  the  result  of  the  vertical  convection  by  the  internal 
waves  of  the  potential  sound-speed  gradient,  i.e.,  the  actual  sound-speed  gradient  minus  the 
adiabatic  gradient 

8c(x,  t )  =  O.c0(z)]p  £(x,  t ),  (35) 

where 

18,  c0(z)]p  =  8,  c0(z)  -  18,  c0(z) ladiabnic  ■  (36) 

The  potential  sound-speed  gradient  is  convected  rather  than  the  measured  gradient  because 

internal  waves  are  oscillations  about  the  adiabatic  background  or  reference  state,  i.e.,  internal 

waves  exist  because  of  entropy  stratification.  With  some  phenomenology  Munk12  and  Munk 

and  Zachariasen3  showed  that 
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Id,  c0(r)]p  -  C0(z)  «J(z)  ,  (37) 

where  g  is  the  acceleration  due  to  gravity,  n  (z)  is,  of  course,  the  bouyancy  frequency  and  /x  is  a 
dimensionless  number  that  can  be  expressed  in  terms  of  certain  numbers  which  describe  frac¬ 
tional  changes  in  several  physical  quantities  with  changes  in  depth.  For  typical  values  of  these 
parameters  they  found  n  —  24.5.  Since  n  is  subject  to  some  variation  this  value  should  only  be 
considered  as  representation.  Desaubies,13  for  example,  found  /x  —  11.6  for  the  sea  water  in 
the  vicinity  of  Cobb  Seamount.  The  fractional  change  in  the  sound  speed  due  to  the  presence 
of  internal  waves  follows  from  Eqs.  (35)  and  (37) 

^  (x,  f)  -  nHz)  £(x,  l).  (38) 

C0  8 

The  mean-square  fractional  sound-speed  fluctuation  can  be  obtained  from  Eqs.  (16)  and  (38) 

-  n0n3(z )  <{o>.  (39) 

r 

Just  as  with  £  we  set 


< 


8  c  ,  x 
—  (x,  /) 
Co 


and  define  an  extrapolated  mean-square  fluctuation  <(8 c/c)q> 


where 


2\l/2 


”o2M 

g 


<Co2>,/2- 


(42) 


For  the  values  n0  “  3  cph,  #x  *  24.5,  and  <{o>l/2  *■  7.3  m  we  obtain 

2\l/2 


5  x  10 


(43) 
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Again,  {  is  a  zero-mean,  Gaussian  random  variable.  Since  the  fractional  sound-speed 
fluctuation  (6c/c0)(x.  /)  is  equal  to  a  non-random,  slowly  varying  function  of  depth  times  £,  it 
also  is  a  zero-mean,  Gaussian  random  variable.  Hence  its  complete  statistical  description  fol¬ 
lows  from  the  correlation  function 


—  (x,  t )  ^  (x\ 

c0  Co 


—  (x,  r)  —  (x'.  n 

Co  C  0 


n  (z) 

7 

6c 

"0 

\ 

c 

This  function  can  be  directly  obtained  from  Eqs.  (31),  (38),  (41)  and  (42) 

.  0 

x  X  J0  dm  W (J) G (w ,n (z))  J0(kjR)  cos[k,(z)(z  -  z’)l  cos[ca(/  -  z')l  .  (44) 

7-1 

Our  goal  here  has  been  the  derivation  of  Eq.  (44).  With  it  one  has  all  that  is  needed  of  the  GM 
model  in  order  to  describe  the  effects  of  internal  waves  on  the  sound-speed  profile. 
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